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摘 要 : 城市 高 速 发 展 引 起 了 区 域 土地 利用 格局 的 改变 ,不 仅 影响 生态 环境 质量 ,而 且 还 对 生态 系 


统 服 务 的 空 


2020 年 净 初 级 生产 力 (Net primary ge N 


的 基础 上 ,将 贝 叶 斯 网 络 和 生态 系统 服务 相 结 


础 上 ,评估 了 4 种 生态 系统 服务 需要 优化 的 区 kA. 为 区 域 经 
明 :(1) 水 文 评价 模型 (Soil and water assessment tool, SWAT) 模 型 能 较 准 确 地 模拟 区 域 的 径流 量 。 


间 格 局 产生 影响 ,在 此 基础 上 优化 生态 系统 服务 显得 至 关 重 要 。 在 泾 河流 域 2000 一 
ue! 、 农 业 生 产 力 .土壤 保持 和 产 水 服务 空间 评估 


合 ,在 关键 变量 子 集 和 和 可视化 的 最 优 状 态 子 集 的 基 
济 和 生态 和 谐 发 展 提 供 参 考 。 结 果 表 


通过 模拟 值 和 观测 值 对 比分 析 ,该 模型 具有 和 较 高 的 决定 性 系数 (R>0.6) 和 纳什 效率 系数 (NSE> 


0.5), 可 为 进 一 


步 评 估 产 水 服务 提供 保障 ,(2) 2000 一 2020 年 泾 河流 域 4 种 生态 系统 服务 的 时 空 差 


异性 较为 显著 。 在 时 间 尺 度 上 ,4 种 生态 系统 服务 均 呈 现 波动 中 上 升 的 趋势 ,在 空间 尺度 上 呈现 较 


为 稳定 的 变化 趋势 (3) 通过 对 4 种 生态 系统 服务 优化 区 域 进行 全 加 分 析 , 发 现 综 
南部 以 及 环 县 的 零星 区 域 。 研 究 结果 对 指导 优化 区 域 进 行 生 态 系统 的 可 持 


在 彭 阳 县 的 中 部 和 


合 优化 区 域 集中 


续 管理 以 及 改善 生态 系统 的 退化 状况 具有 重要 意义 。 


关 键 词 : 生态 系统 服务 
文章 编号 : 


生态 系统 服务 是 指 人 们 从 生态 系统 中 获得 的 
所 有 受益 ,包括 直接 和 间接 受益 "“”。 由 于 生态 系 
统 服 务 与 人 类 的 生命 和 生存 密切 相关 ,因此 研究 生 
态 系统 服务 对 于 人 类 的 生存 和 地 球 发 展 都 非常 重 
要 。 随 着 科学 技术 水 平 的 提升 , 越 来 越 多 的 学 者 
对 生态 系统 服务 研究 集中 于 供需 关系 、 空 间 权 衡 与 
优化 .生态 福 社 、 服 务 流 等 方面 "”。 其 中 生态 系统 
服务 空间 格局 优化 受到 越 来 越 多 的 学 者 关注 ,学 者 
们 开始 尝试 运用 不 同 模 型 与 方法 进行 深入 探究 ,但 
是 由 于 研究 时 限 不 长 ,所 以 促使 该 方面 研究 处 于 萌 
芽 时 期 。 目 前 多 数 空间 优化 研究 是 以 土地 利用 
空间 格局 为 支撑 ,其 目的 是 实现 一 定 规模 的 需求 和 
规划 的 目标 。 刘 慧 灵 等 ”通过 单位 面积 上 不 同 
土地 利用 的 生态 系统 服务 价值 来 作为 变量 系数 , 建 
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立 生态 效益 函数 ,进而 对 土地 利用 进行 空间 优化 ; 
王 观 清 等 ”向 入 靶 等 "通过 土壤 有 机 碳 的 最 大 储 
量 为 目标 ,确定 生态 效益 函数 ,最 终 达 到 优化 土地 
利用 的 目的 。 李 秀 霞 等 ”通过 设置 经 济 和 生态 的 
双重 目标 ,利用 系统 动力 学 模型 和 多 目标 优化 模 
型 ,对 吉林 省 的 土地 利用 进行 仿真 模拟 和 优化 。 上 
述 研究 大 多 数 是 对 土地 利用 进行 空间 优化 ,然后 在 
以 优化 的 土地 利用 为 基础 ,进行 生态 系统 服务 的 空 
间 制 图 和 分 析 ,不 仅 忽 略 影响 生态 系统 服务 的 其 他 
影响 因素 ,而且 没 有 真正 意义 上 实现 生态 系统 服务 
的 空间 优化 , 即 在 研究 区 范围 内 ,如何 合 理 和 优化 
配置 自然 因素 和 土地 利用 等 因素 ,使 该 区 域 的 生态 
系统 服务 能 力 最 强 ,如 何 选取 不 同 生 态 系 统 服务 的 
影响 因素 以 及 量化 关键 因子 对 生态 系统 服务 的 影 
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响 程度 ,从 而 达到 优化 区 域 生 态 系统 服务 的 目的 ， 
正 是 该 研究 所 解决 的 问题 。 

贝 叶 斯 网 络 是 一 种 大 规模 的 数据 分 析 工 具 , 具 
有 强大 的 逻辑 分 析 能 力 。 该 网 络 采 用 图 形 化 建 模 
的 方法 ,将 变量 之 间 的 相互 作用 以 网 络 的 形式 表达 
出 来 ,并 构建 一 个 系统 概率 模型 ,将 贝 叶 斯 网 络 
和 生态 系统 服务 相 结 合 , 可 以 很 好 地 为 生态 系统 服 
务 的 优化 管理 提供 决策 文 持 ”。Fox 等 研究 了 贝 
叶 斯 网 络 在 水 质 和 自然 资源 管理 中 的 价值 ;Landuyt 
等 “利用 贝 叶 斯 网 络 评 估 池 塘 综 合 管理 的 机 会 和 
风险 ;Dang 等 ”将 生态 系统 服务 和 自然 环境 因子 纳 
和信 贝 叶 斯 网 络 模 型 ,通过 不 同 指标 设置 不 同 的 研究 
情景 ,为 农业 发 展 提供 概率 支持 。 综 上 看 出 贝 叶 斯 
网 络 可 以 很 好 地 和 生态 系统 服务 相 结 合 , 从 而 来 实 
现 研 究 目标 。 

泾 河流 域 位 于 黄土 高 原 中 部 ,在 六 盘山 和 子午 
岭 之 间 , 对 于 此 区 域 的 研究 大 多 集中 在 生态 和 自然 
环境 方面 。Chen 等 “使 用 主 成 分 分 析 和 相关 分 析 
来 分 析 泾 河流 域 的 归 一 化 植被 指数 (Normalized dif- 
ference vegetation index, NDVI) 和 降水 量 , 发 现 植被 
的 动态 变化 与 气候 变化 、 洪 水 和 人 类 活动 密切 相 
关 ;Xie 等 ”分 析 了 降水 -景观 -径流 变化 和 相互 作 
用 ;Guo 等 ”使 用 滑动 偏 相 关系 数 法 分 析 降 十 和 径 
流 关系 的 变化 ,并 使 用 双 累 积 曲 线 法 验证 其 分 析 ; 
HIRE O 通过 参与 式 社区 评估 法 对 泾 河流 域 景观 
管理 中 存在 的 相关 问题 进行 次 入 衣 析 。 从 现 有 研 
究 中 看 出 对 于 泾 河流 域 生 态 系统 服务 的 研究 还 处 
在 初级 阶段 ,并 且 随 着 社会 经 济 的 快速 发 展 ,土地 
利用 格局 发 生 很 显著 的 变化 ,最 终 也 会 导致 生态 系 
统 服务 的 格局 发 生变 化 ,如 何 通 过 现 有 的 土地 利用 
格局 来 优化 生态 系统 服务 也 显得 至 关 重 要 ,因此 ， 
深入 研究 泾 河流 域 生态 系统 服务 以 及 空间 格局 优 
化 会 对 区 域 的 生态 -经 济 和 谐 发 展 黄 定 基础 。 本 文 
在 生态 系统 服务 空间 评价 的 基础 上 ,引入 贝 叶 斯 网 
络 模型 ,通过 筛选 关键 因子 和 最 优 状态 子 集 的 方 
式 , 来 获取 4 种 生态 系统 服务 的 空间 优化 区 域 , 为 区 
域 生态 健康 持续 发 展 葛 定 基础 。 


1 数据 与 方法 


1.1 研究 区 概况 
泾 河 发 源 于 泾 源 县 弃 源 乡 , 源 区 地 处 六 盘山 东 


苞 , 自 西北 向 东南 流 经 固原 市 平凉 市 、 庆 阳 市 . 咸 
阳 市 等 38 个 县 市 ,于 咸阳 市 高 陵 县 陈 家 滩 村 汇 人 
渭河 。 泾 河流 域 地 理 位 置 介 于 106?14'~108?"42" 工 、 
34*46' -37?19' N Z [E] ,地 跨 陕 西 省 .甘肃 省 .宁夏 回 
族 自治 区 (图 1) ,流域 面积 达 45421 kms KREE 
平均 气温 8 % ,年 降水 量 350~650 mm, 降水 主要 集 
中 在 夏季 。 泾 河流 域 地 势 较为 平坦 ,耕地 面积 约 占 
1/3 ,农业 发 达 ; 矿 藏 量 较为 丰富 ,经济 开 发 潜力 大 。 
泾 河流 域内 黄土 高 原 区 域 和 泾 囊 渠 灌区 域 ,是 陕西 
省 主要 商品 粮油 基地 。 除 此 之 外 ,还 有 煤炭 石油 、 
天 然 气 建筑 材料 等 矿产 ,如 长 庆 油田 。 但 长 期 过 
度 开发 利用 和 气候 变迁 ,导致 流域 内 生态 环境 十 分 
脆弱 ,其 中 水 资源 缺乏 已 成 为 研究 区 社会 经 济 以 及 
生态 发 展 的 主要 制约 因素 。 
1.2 数据 来 源 

本 文 主要 数据 包括 遥感 数据 (土地 利用 数据 ) 
和 非 遥 感 数据 (地 理 基 础 数据 .高 程 数 据 .气象 数 
据 .土壤 类 型 数据 和 水 文 数据 等 ) ,具体 信息 见 表 1 
所 示 。 
1.3 生态 系统 服务 模型 
1.3.1 农业 生产 力 服务 模拟 ”首先 ,从 2000 一 2020 
年 地 表 和 覆盖 数据 中 获取 耕地 资源 数据 作为 计算 区 


图 例 

A 张 家 山 水 文 站 ”海拔 高 度 /m 

。 县 级 行政 中 心 m 2908 
219 


一 泾 河 支 流 
一 泾 河 干流 
图 1 研究 区 概况 
Fig. 1 Overview of the study area 
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表 1 数据 来 源 


Tab. 1 Data sources 


数据 类 型 数据 描述 数据 来 源 

基础 地 理 信 息 数据 行政 区 划 河流 .道路 等 国家 基础 地 理 信 息 中 心 (http:Wwww.ngcc.cn/ngcc/) 

高 程 数 据 来 源 于 地 理 空间 数据 云 ,空间 分 辨 率 为 30 m 地 理 空间 数据 云 (http://www.gscloud.cn/) 

土壤 数据 第 二 次 全 国土 地 调查 南京 土壤 研究 所 提供 的 1:100000 土 ”国家 青藏 高 原 科学 数据 中 心 (https://data.tpde.ac.cn/zh- 
壤 数 据 ,空间 分 辨 率 为 1 kmx1 km hans/) 

土地 利用 数据 来 源 于 GlobalLand 30, 全 球 地 理 信 息 公 共产 品 , 空 间 分 辨 全球 地 理 信息 公共 产品 (http://www.globallandcover. 
率 为 30 m com/) 

气象 数据 2000 一 2020 年 逐日 数据 ,指标 包括 降水 量 .最 高 气温 .最 。” 国家 气象 科学 数据 中 心 (http://data.cma.cn/) 
低 气温 .风速 .相对 湿度 数据 

水 文 数据 2000 一 2020 年 张 家 山水 文 站 逐日 径流 数据 张 家 山 水 文 站 (http://61.163.88.227:8006/hwsq.aspx) 

统计 数据 2000 一 2020 年 人 口 密度 以 及 GDP 等 统计 数据 统计 年 鉴 


域 农 业 生 产 潜力 的 基础 ;其 次 是 基于 气象 站 点 数 
据 , 通 过 ANUSPLIN 软件 获取 2000 一 2020 平 均 气 温 
和 年 降水 量 数据 ;最 后 ,采用 桑 斯 维特 纪念 模型 ” ， 
计算 区 域 2000 一 2020 年 农业 生产 力 ,公式 如 下 : 

CP-PvXI,, (1) 
式 中 :CP 为 以 耕地 为 基础 的 农作物 生产 力 (t:hm”); 
Pv 为 以 年 平均 蒸 散 量 为 自 变量 计算 出 的 植物 干 物 
质 总 量 (g.m2.a ;为 土地 利用 水 平 , 是 由 中 国 农 
用 地 分 级 规定 决定 的 320。Pv 的 计算 公式 如 下 : 


—0.0009695 x (v — 20) 


Py 230000 xl -e ] (2) 
1.05, 
s VS I. (3) 
[1 +(.05pily 
12300 + 251 + 0.05? (4) 


式 中 :2 为 年 平均 实际 蒸 散 量 (mm ) jp 为 年 降水 量 
(mm) ;为 年 平均 最 大 薰 散 量 (mmy) ;为 年 平均 气温 
(OS 
1.3.0 产 水 服务 模拟 ”水文 评价 模型 (Soil and wa- 
ter assessment tool, SWAT) 可 以 对 流域 内 一 系列 复 
杂 的 物理 过 程 进行 模拟 ,模型 中 的 水 量 平衡 始终 是 
流域 内 所 有 过 程 的 驱动 力 他 。 并 且 此 模型 可 以 基 
于 水 量 平衡 模拟 每 个 水 文 啊 应 单元 的 地 表 径 流量 
和 洪峰 流量 ,水 量 平衡 方程 如 下 

SW, =SW,+ X (Ra, -0 E, -Won -0,) (5) 

i=l 


式 中 :SW, 为 最 终 土 壤 含水 量 (mm ) 六 为 时 间 天 数 ; 
SW 为 第 i 天 的 初始 土壤 含水 量 (mm) ;Ri 为 第 i 天 
的 降水 量 (mmy) ;0O. 必 为 第 ;天 的 地 表 径 流量 (mm ) ; E, 
HF i RZS (mm) ;本 .为 第 ;天 通过 土壤 剖面 
进入 包 气 带 的 水 量 (mm) ;0 为 第 i 天 的 回归 流量 


(mm). 


1.3.3. 净 初 级 生产 力 (NPP) 服 务 模 拟 ”参考 朱文 泉 
等 中 的 遥感 估算 模型 (Carnegie-ames- stanford ap- 
proach, CASA ) 计 算 植被 NPP, 由 植物 的 光合 有 效 辐 
射 (Absorbed photosynthetic active radiation, APAR ) 和 
实际 光 能 利用 率 ( e )2 个 因子 表示 ”。 计 算 公 式 如 
Tes 

NPP(x,t)= APAR(x,t) X e(x,t) (6) 
XX'P NPP(w.t) 9 IC x dE t H N ATIR S ERRO 
级 生产 力 (g Cm?)  APARGOD JU x TE t H WHS 
光合 有 效 辐射 (MJm ”>) ; e(t ARC x TE t H KER 
光 能 利用 率 (gC.MJ )。 
1.3.4 土壤 保持 服务 模拟 ， 采 用 通用 土壤 流失 方程 
估算 土壤 保持 量 ,20 志 纪 50 年 代 由 Wischmeier $E 
提出 ,基于 理论 以 及 大 量 实地 观测 数据 统计 分 析 的 
基础 上 ,对 各 因子 都 有 具体 解释 。 通 用 土壤 流失 方 
程 的 表达 如 下 : 

A-RXKXLXSX(1-CxP) (7) 
式 中 :4 为 土壤 保持 量 (t.hm”) ;RR 为 降雨 侵蚀 力 因 子 
(MJ: mm-hm^-*h^*a ;天 为 土壤 可 亿 性 因子 (tn 
MJ mm ) 记 和 3S 分 别 为 坡 长 和 坡度 因子 ;C 为 地 表 
植被 覆盖 及 经 营 管理 因子 ;为 水 土 保持 措施 因子 。 
1.4 基于 BBNs 的 农田 生态 系统 服务 空间 格局 优化 
14.1 贝 叶 斯 网 络 模型 ， 贝 叶 斯 网 络 又 称 信 度 网 
络 , 是 Bayes 方 法 的 扩展 ,目前 不 确定 知识 表达 和 推 
理 领域 最 有 效 的 理论 模型 之 一 ”。1988 年 Pearl 提 
出 后 ,已 经 成 为 近 几 年 来 研究 的 热点 。 一 个 贝 叶 斯 
网 络 是 一 个 有 向 无 环 图 (DAG), 由 代表 变量 节点 及 
连接 这 些 节 点 有 向 边 构 成 ”。 

其 中 条 件 概 率 表 (CPT) 表 达 了 父 节 点 X 与 子 节 

点 了 之 间 的 关系 强度 , 表 中 每 一 行 由 父 节 点 的 状态 


202208.00076v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


余 玉 洋 等 : 基于 贝 叶 斯 网 络 的 生态 系统 服务 空间 格局 优化 一 一 以 泾 河流 域 为 例 


组 合 和 条 件 概 率 POlw) 组 成 , 即 父 节 点 X 在 离散 状态 
x 下 , 子 节 点 对 应 状态 y 发 生 的 概率 。 子 节点 的 概率 
分 布 由 父 节 点 的 概率 分 布 及 条 件 概率 表决 定 , 对 于 
没有 父 节 点 的 节点 ,其 概率 分 布 为 先 验 概率 分 布 P 
(加 。 所 有 节点 的 条 件 概率 分 布 相 乘 得 到 贝 叶 斯 网 
络 中 所 有 变量 的 联合 概率 [ 式 (8)]。 贝 叶 斯 网 络 能 
够 有 效 计算 网 络 中 任意 节点 概率 分 布 的 基础 ,为 本 
文 优 化 生态 系统 服务 格局 提供 了 理论 "1。 
P(X,.X,, X) - ET. ,P[Xdparents(X;] (8) 
XP P(X, X, XAA AARE ES ES A) B L3 
种 分 布 的 随机 变量 具有 1,2, n T ESRB in 2 BG 
PEERK; (XX, ,XX,) 为 随机 变量 。 
以 4 种 生态 系统 服务 量化 模型 为 依据 ,选择 与 
其 相关 的 影响 因子 作为 节点 ,构建 生态 系统 服务 的 
贝 叶 斯 概念 网 络 (图 2)。 借 助 AreGIS 软件 ,离散 市 


点 对 应 的 栅 格 图 层 , 根 据 变量 的 实际 分 布 情况 , 参 
考 学 者 对 贝 叶 斯 网 络 栅 格 分 层 的 方法 ”2 ,运用 Arc- 
GIS 软件 的 自然 间断 法 将 上 述 每 个 变量 栅 格 数据 分 
为 4 个 等 级 :最 高 .高 .中 、 低 ,具体 离散 化 分 级 标准 
如 表 2 所 示 。 提 取 离 散 数据 值 到 属性 表 并 将 概率 条 
件 输入 贝 叶 斯 网 络 模型 , 便 可 得 到 泾 河 流域 生态 系 
统 服 务 的 贝 叶 斯 网 络 。 其 他 节点 的 条 件 概 率 表 可 
由 条 件 概 率 公 式 [ 式 (9)] 计 算 获得 。 

P(BIA)= P(AB)/P(A) (9) 
式 中 :P(AB) 为 事件 A 与 B 的 联合 概率 , 即 2 个 事件 
共同 发 生 的 概率 ;P(A) 为 A 的 先 验 概率 , 即 事件 A 发 
生 的 概率 ; 同 理 ,P(B) 为 事件 B 发 生 的 概率 。 
1.4.2 农田 生态 系统 服务 格局 优化 ”构建 出 泾 河流 
域 生态 系统 服务 贝 叶 斯 网 络 模型 ,计算 出 条 件 概率 
表 和 敏感 性 因子 ,确定 出 关键 变量 关键 状态 子 集 。 


(a) 数据 图 层 (b) 属性 表 
k i A 像 元 点 降水 量 ”植被 ”坡度 RE tA ”太阳 辐射 气温 /PC 
图 层 2 = ME py > “从 [mm 覆盖 度 (O) Emm 用 水 平分 级 /MJ-m? 
Rs NEU m 1 55840 >0.46 >7.0 >5362.0 3 255760  »90 
图 层 n 图 层 n Lal 
T, 状态 分 级 
(c) 贝 叶 斯 模型 ES , CERES 
耕地 42.4 L3 因素 概率 /% 
林地 11.0 salt 最 高 54.5 
草地 ”44.0 最 高 0.5 高 42.9 
水 域 ”0.5 高 18.4 中 03 
建设 用 地 2.1 ie E 低 23 
MEE 0.0 IARE 
on 最 高 35 
最 高 16.5 M a 
高 9.7 i 
中 32.9 低 62.0 im 
低 409 概率 /% 
=n 太阳 辐射 最 高 23.3 
气温 概率 /% | 降水 量 概率 /24| 概率 /% 高 310 
最 高 12.7 最 高 最 高 313 中 30.7 
中 51.5 低 "n 中 13.1 
低 0.0 5 低 45 
(d) 条 件 概 率 表 4 
土壤 保持 量 概 率 /% 
一 植被 覆 蒸 散 ” 土地 利用 “太阳 
像 元 点 降水 量 cogn BUE 发 量 水 平分 级 mi T 最 高 高 中 低 


1 最 高 最 高 高 高 


高 高 高 3.5 8.9 25.6 62.0 


Fig.2 Schematic diagram of Bayesian belief 


注 :NPP 为 净 初 级 生产 力 。 下 同 。 
图 2 生态 系统 服务 贝 叶 斯 网 络 构建 原理 


network construction for farmland ecosystem service 
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FHEWY 


GE: NPP 为 净 初 级 4 


表 2 因子 状态 分 级 


Tab. 2 State classification of factors 


因子 


蒸 散发 量 /mm 


土地 利用 水 平分 级 


k 


É Hr 


植被 覆 


B 
X 


水 土 保持 措施 因素 


土壤 类 型 


太阳 辐射 /MJ.m” 


降水 量 /mm 


Ji REI) 


气温 /C 


土地 利用 类 型 


土壤 保持 /thm? 


AALS? Ft hm? 


NPP/g C* m? 


产 水 量 /mm 


E 产 力 。 


状态 名 称 


kid kid kid a in Ek 


up 


X -EOza ZR -R zm mL xU GRozn mm mm -RE orn Hn om» -m zn mr ox: Gora mm mm Roin Hn om» cm zm AHH 


未 利用 地 


uj 


和 zb + Sm + 


x 


ta 
hi 


实际 值 范围 
2256~3903 
3903~5362 
5362~9827 
9827~65534 
1 
2 
3 
4 
-0.17~0.14 
0.14~0.30 
0.30~0.46 
0.46~0.70 
0.00 
0.00~0.10 
0.10~0.35 
0.35~1.00 
1~4 
4-9 
9-14 
14-20 
4989-5331 
5331-5576 
5576-5732 
5732-5911 
323-422 
422-506 
506-584 
584-693 
0.00-3.50 
3.50~7.00 
7.00~10.50 
<8 
8~9 
9~10 
10~12 


0.00-41.30 
41.30-93.60 
93.60~174.14 
174.14~555.09 
28.64~57.75 
57.75~60.00 
60.00~90.00 
90.00~ 120.00 
0.23~11.97 
11.97~135.89 
135.89~204.47 
204.47~388.90 
0.00~408.61 
408.61~1001.75 
1001.75~2095.70 
2095.70~3361.15 


以 ArcGIS 软件 为 工具 ,可 视 化 子 集 分 布 空间 ,最 终 
绘制 不 同等 级 所 对 应 的 关键 变量 关键 状态 子 集 同 
时 发 生 的 区 域 。 筛 选 关键 变量 关键 状态 子 集 ,包括 
以 下 步 又: 

(1) 确定 各 变量 的 关键 状态 

根据 节点 各 状态 的 概率 及 两 两 节点 状态 组 合 
的 联合 概率 ,计算 每 个 影响 因子 与 各 服务 间 的 条 件 
概率 ,选取 不 同等 级 对 应 各 个 影响 因子 的 条 件 概 率 
最 大 状态 作为 变量 的 关键 状态 。 当 服务 处 在 最 高 、 
高 .中 、 低 4 种 不 同 级 别 时 ,将 对 应 4 种 变量 状态 

A 

(2) 确定 关键 变量 

为 了 评价 贝 叶 斯 网 络 节 点 的 相对 重要 性 ,Neti- 
ca 提供 了 敏感 性 分 析 ”” , 主要 是 用 于 评估 生态 系 
统 服 务 节 点 是 否 会 敏感 地 感知 其 他 影响 因子 节点 
的 变化 ””™。 据 此 确定 生态 系统 服务 关键 变量 关键 
状态 子 集 ,然后 选取 变量 关键 状态 子 集 的 面积 与 研 
究 区 总 面积 比值 最 大 的 情景 作为 生态 系统 服务 最 
优 的 情景 ,将 该 情景 下 的 关键 变量 关键 状态 子 集 选 
定 为 最 优 状 态 子 集 , 然 后 结合 县 区 的 发 展 和 管理 等 
划分 出 生态 系统 服务 的 优化 区 。 


2 结果 与 分 析 


2.1 SWAT 模 型 率 定 与 验证 

SWAT 模 型 有 26 个 与 径流 模拟 相关 的 敏感 性 
子 , 不 同 因 子 影响 各 异 。 本 文 根 据 SWAT-CUP 软 
件 进 行 模拟 率 定 ,通过 了 检验 确定 水 文 参 数 的 敏感 
性 相对 显著 性 ,P 值 确定 其 显著 性 特征 人 。 进 而 
选取 相关 参数 ( 表 3) ,通过 调整 参数 ,误差 能 够 控制 
在 一 定 范围 之 内 。 

由 于 张 家 山 水 文 站 观测 数据 限制 ,本 文 以 
2000 一 2010 年 实测 数据 和 模拟 数据 进行 率 定 ,根据 
不 断 调 整 得 到 的 参数 ,获得 模拟 效果 最 佳 的 一 组 数 
据 。 如 图 3a 所 示 ,其 中 决定 性 系数 ( 尼 ) 为 0.74, 纳 什 
效率 系数 (NSE ) 为 0.73 ,都 大 于 0.5, 表 示 模 拟 值 与 
实测 值 具有 较 高 的 相关 性 ,模型 的 模拟 结果 是 可 靠 
的 。 模 拟 2011 一 2020 年 泾 河流 域 的 月 径流 量 ,与 张 
家 山水 文 站 数据 进行 比 对 进行 验证 ,可 以 得 出 尼 为 
0.78( 图 3b) ,径流 实测 与 模拟 值 具有 较 高 的 相关 性 ， 
模拟 结果 是 可 信 的 , 据 此 模拟 2000 一 2020 年 泾 河流 
域 产 水 服务 。 
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表 3 参数 敏感 性 分 析 


Tab.3 Sensitivity analysis result of parameters 


参数 名 称 相对 显著 性 (7) 显著 性 检验 (P) 参数 范围 最 优 校准 值 
土壤 饱和 导 水 率 2.5450 0.0112 -0.9~0.9 -0.8425 
降雪 温度 -2.2385 0.0256 -20-20 7.9700 
主 河道 水 力 传导 度 1.8577 0.0638 0~150 58.8517 
地 下 水 清 留 系数 -1.7717 0.0771 0~500 390.1691 
EAE AME RA -1.2422 0.2148 0.01~1.00 0.0294 
土壤 有 效 含 水 率 -1.1148 0.2655 -0.5~0.5 0.3959 
土壤 湿 密度 1.0494 0.2945 -0.5~0.5 -0.4496 
基 流 系数 0.5965 0.5512 0-1 0.9276 
地 下 水 再 蒸发 系数 -0.5050 0.6138 0.02~0.20 0.1400 
浅 层 地 下 水 径流 系数 0.0329 0.9737 0~5000 2409.2715 

250 250 


(a) 率 定 期 一 实测 值 一 模拟 值 


Pp 200 R1-0.74 
ES 150 NSE-0.73 
xe 100 
V so 
We Tu 9 tl 
0 P u . a n > 
CO SS HDA AD 
X ALA ^ BY QS a ^ 
SS SS S s KM SS S SS SS SS gy 
"m e fy" aun v C^ y yt fy yc Cy cv) 
日 期 /年 -月 ) 


(b) 验证 期 


一 实测 值 一 模拟 值 
R=0.78 


150 NSE-0.78 
100 
WV so 
0 
d d od & dx doa od d 
S. S. S as. Ss Ka ooo 
日 期 (年 -月 ) 


注 :RR 为 决定 性 系数 ;NSE 为 纳什 效率 系数 。 
图 3 率 定 期 和 验证 期 逐 月 径流 量 模拟 值 与 实测 值 


Fig. 3. Simulated and measured values of monthly runoff during the regular period and the verification period 


2.2 生态 系统 服务 时 空格 局 分 析 

生态 系统 服务 功能 之 间 相 互 作用 与 影响 ,对 流 
域 生态 格局 发 展 与 优化 具有 显著 影响 。 本 研究 通 
过 模型 模拟 了 2000—2020 年 泾 河 流域 产 水量 、 
NPP 农业 生产 力 和 土壤 保持 4 种 生态 系统 服务 量 
的 空间 分 布 (图 4)。 区 域 产 水 量 在 时 间 上 呈 增 长 趋 
势 ,2000 年 平均 产 水 量 为 489.6 m, 而 2020 年 平均 
产 水 量 为 729.6 m。 在 空间 上 呈现 北部 和 西南 部 
高 ,中 间 低 的 特征 ,高 值 区 域 主要 集中 在 环 县 . 千 阳 
县 和 凤翔 县 等 区 域 。 此 外 , 低 值 区 主要 集中 在 西部 
的 泾 源 县 以 及 北部 的 盐池 县 。 但 是 2020 年 产 水 量 
高 值 区 域 较为 分 散 , 主 要 呈现 中 间 低 ,四 周 高 的 分 
布 格局 ,主要 是 因为 2020 年 相 较 于 其 他 年 份 降水 重 
心 有 所 转移 ,并 且 受 季风 天 气 的 影响 ,导致 区 域 降 
水 量 年 际 差 异 明显 ,所 以 在 一 定 程度 上 影响 产 水 量 
的 空间 分 布 。NPP 在 年 际 变化 中 呈现 增长 的 趋势 ， 
2000 年 NPP 的 平均 值 为 185.1 g Cm”, 总 量 为 7.7x 
10°t, mi 2020 fE NPP 的 平均 值 为 461.9 g C* m^, iti 


为 2.4x10't, 总 量 相 较 于 2000 年 大 约 增加 了 3 倍 。 
NPP 在 空间 上 呈现 西北 向 东南 逐渐 递减 的 趋势 ,高 
值 区 域 集 中 在 流域 两 侧 子 午 岭 山系 及 西部 六 盘山 
地 区 。 低 值 区 主要 原因 是 植被 以 典型 草原 为 主 , 并 
且 逐 渐 东 南方 向 退化 ,沙漠 化 问题 较 严 重 。 农 业 生 
产 力 在 时 间 上 从 由 2000 年 的 280.9x10"t, 增 长 到 至 
2020 年 的 368.3x107t; 在 空间 上 自 西向 东 呈 增加 趋 
势 ,南北 向 旦 "马鞍 形 ” 分 布 , 空 间 分 布 差异 较 小 。 
高 值 分 布 面积 有 所 增加 ,集中 分 布 在 陕西 省 长 武 
县 \ 淳 化 县 以 及 宁夏 回族 自治 区 泾 源 县 。 土 壤 保 持 
从 时 间 尺 度 上 看 ,总 体 呈 波动 变化 趋势 ,但 年 际 差 
异 较 大 ,土壤 保持 从 2000 年 的 3.0x10t 增 加 到 2020 
4E. 5.3x 10 t, 土壤 保持 呈现 显著 增长 趋势 。 从 空 
间 分 布 特征 看 出 ,土壤 保持 强度 空间 分 布 存 在 一 定 
差异 ,整体 呈现 出 北部 .中 部 低 ,西南 部 和 东南 部 高 
的 分 布 格局 。 土 壤 保持 强度 低 值 区 主要 分 布 在 流 
域 北 部 和 中 部 ,北部 为 地 形 起 伏 变 化 的 黄土 丘陵 
区 ,加 之 年 降水 量 低 于 区 域 平均 值 ,植被 生长 水 热 
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(a) 2000 年 产 水 量 (b) 2005 年 产 水 量 
盐池 县 局 


(c) 2010 年 产 水 量 


(d) 2015 年 产 水量 


(e) 2020 年 产 水 量 


3361.2 


图 例 
c 县 域 边界 
NPP/g C-m? 


(k) 2000 年 农业 生产 力 ”(1) 2005 年 农业 生产 力 ”(m) 2010 年 农业 生产 力 ”(n) 2015 年 农业 生产 力 (0) 2020 年 农业 生产 力 


(P) 2000 年 土壤 保持 (q) 2005 年 土壤 保持 


(r) 2010 年 土壤 保持 


(s) 2015 年 土壤 保持 


(t) 2020 年 土壤 保持 


图 例 
口 县 域 边界 
土壤 保持 /thm? 


图 4 弃 河 流域 生态 系统 服务 量 空间 变化 与 分 布 


Fig.4 Change and distribution of ecosystem Services in Jing River Basin 


条 件 很 差 ,导致 该 地 区 植被 覆盖 度 较 低 。 高 值 区 主 
要 分 布 在 流域 西南 部 和 东南 部 ,部 分 区 域 大 面积 ， 
林地 ,海拔 较 高 ,并 且 规划 了 和 森林 保护 区 ,人 类 活动 
对 地 表 的 干扰 小 。 
2.3 生态 系统 服务 空间 格局 优化 

贝 叶 斯 网 络 模型 以 概率 相关 理论 为 基础 ,得 到 
研究 区 不 同 生态 系统 服务 关键 变量 最 优 状 态 子 集 ， 
进而 筛选 出 需要 优化 的 区 域 。 对 于 泾 河流 域 生 态 
系统 来 说 ,人 为 因素 和 自然 环境 的 双重 作用 ,对 生 
态 系统 服务 的 影响 不 容 小 舰 。 不 同 的 生态 系统 服 


FB ,筛选 不 同 的 关键 变量 ,得 到 不 同 生态 系统 服务 
的 条 件 概 率 表 , 如 表 4~7 所 示 。 根 据 表 4 所 示 , 土 地 
利用 水 平和 蒸 散发 是 影响 农业 生产 的 关键 变量 , 当 
农业 生产 力 出 现 最 高 (0.353) 、 高 (0.422) 、 中 等 
(0.620) 和 低 水 平 (0.491) 概 率 时 ,所 对 应 的 土地 利 
用 水 平和 蒸 散 发 量 分 别处 于 不 同 程度 的 水 平 。 对 
比 得 到 土地 利用 水 平 处 在 高 级 水 平 即 农业 用 地 水 
平 , 蒸 散发 量 处 于 高 值 区 间 5362-9827 mm 时 ,农业 
生产 力 达 到 最 大 可 能 概率 (0.620) ,处 于 中 值 范 围 
(57.7~60.0 t- hm?) ,对 其 进行 空间 可 视 化 展示 (图 
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表 4 农业 生产 力 节点 条 件 概率 
Tab.4 Conditional probability of the node crop 


用 ,能 在 一 定 程度 上 增加 地 表 径 流量 ,最 终 提 升 区 
域 的 产 水 服务 能 力 。 当 土地 利用 类 型 为 未 利用 地 ， 


ee 太阳 辐射 .气温 和 植被 覆盖 度 都 为 最 高 水 平时 ,NPP 
土地 利用 水 平 ” 菩 散 发 量 —— pe 服务 达到 最 大 可 能 概率 (0.878) ,处 于 中 值 范围 
we (11.97~135.89 g C ^m?) ,优化 区 域 主要 集中 在 彭 阳 
最 高 最 高 0.353 0.002 0284 0.361 : e 
s nj 县 和 华 池 县 ,镇 原 县 、 环 县 中 部 和 南部 ,吴起 县 南部 
I] z 0.220 0.032 0.620 0.119 i 、 m 
P A add dtii -0299 quid) 和 合 水 县 北部 ,提升 区 域 的 植被 覆盖 率 对 提高 NPP 
低 中 0.108 0.422 0.282 0.188 FRA A PRE AY FM, 38 zd SR EES PAY TSE ,影响 土 


5) ,优化 区 域 主要 分 布 在 环 县 WA BAR 
阳 县 以 及 镇 原 县 等 北部 黄土 区 ,对 这 些 地 区 加 大 农 
业 生 产 的 投资 和 管理 力度 ,能 够 较为 迅速 地 提高 泾 
河流 域 整体 的 农业 生产 力 。 当 土地 利用 类 型 为 草 
地 ,土地 利用 水 平 处 在 中 级 水 平 , 降 水 量 处 于 最 高 
值 584~693 mm 时 ,气温 小 于 8 % 时 , 产 水 量 达 到 最 
大 可 能 概率 (0.805) ,处 于 高 值 范围 (408.61~ 
1001.75 m^) ,优化 区 域 主要 集中 在 环 县 的 中 南部 以 
及 彭 阳 县 和 镇 原 县 附近 ,对 区 域 水 资源 进行 高 效 利 


二 保持 的 关键 因素 是 降水 量 .坡度 ,水土 保持 措施 
因素 .植被 和 土壤 类 型 ,降水 量 概率 是 最 重要 的 
Ro MERKA .降水 量 坡度 以 及 水 土 保持 措施 
因素 位 于 中 级 水 平 ,植被 覆盖 度 为 最 高 水 平时 ,所 
对 应 的 土壤 保持 出 现 最 高 水 平 的 概率 最 高 ,达到 
0.824 ,该 子 集 主 要 分 布 在 欧阳 县 和 华 池 上 县. 庆 城 县 
北部 和 镇 原 县 以 及 环 县 南部 ,对 该 区 域 增加 植被 覆 
盖 , 抚育 和 保护 森林 资源 ,加 强生 态 环境 的 建设 和 
投入 ,有 效 控制 水 土 流失 带 来 的 损失 ,会 对 区 域 的 
水 土 保持 服务 起 到 正 向 的 作用 。 


表 5 产 水 节点 条 件 概率 
Tab.S Conditional probability of the node water yield 


土地 利用 水 平 降水 量 气温 土地 利用 类 型 — IEEE 
最 高 高 中 低 

最 高 最 高 高 未 利用 地 0.159 0.020 0.205 0.616 

中 最 高 低 未 利用 地 0.065 0.138 0.579 0.217 

中 最 高 低 草地 0.031 0.805 0.084 0.080 

最 高 最 高 最 高 草地 0.488 0.124 0.191 0.196 

#6 NPP 节 点 条 件 概率 
Tab.6 Conditional probability ofthe node NPP 
降水 量 气温 土地 利用 类 型 太阳 辆 射 植被 覆盖 度 TE E 
最 高 高 中 低 

中 最 高 未 利用 地 最 高 最 高 0.060 0.024 0.878 0.039 

低 低 未 利用 地 低 低 0.072 0.043 0.438 0.446 

中 最 高 草地 最 高 最 高 0.149 0.469 0.019 0.363 

最 高 最 高 耕地 最 高 最 高 0.571 0.154 0.160 0.115 

表 7 土壤 保持 节点 条 件 概率 
Tab. 7 Conditional probability of the node soil conservation 
土壤 类 型 降水 量 植被 覆盖 度 坡度 土壤 保持 措施 因素 = = a 

最 高 高 中 低 
最 高 最 高 中 低 最 高 0.126 0.056 0.320 0.499 
最 高 最 高 最 高 低 最 高 0.191 0.135 0.458 0.215 
最 高 最 高 最 高 中 高 0.108 0.496 0.035 0.362 
中 中 最 高 中 中 0.824 0.007 0.139 0.029 
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最 后 对 4 种 生态 系统 服务 的 优化 区 域 进 行 空间 
LAT ,得 到 综合 的 优化 像 元 (图 5)。 优 化 区 域 
主要 集中 在 彭 阳 县 的 中 部 和 西南 部 以 及 环 县 的 零 
星 区 域 , 部 分 区 域 应 该 采取 相应 的 生态 环境 保护 措 
施 ,加快 生态 .民生 和 多 功能 林业 的 发 展 ,对 于 草原 
生态 要 坚持 保护 与 建设 并 举 ,保护 优先 的 原则 , 设 
置 一 定 的 惩 奖 措施 ,在 一 定 程度 上 提升 居民 对 于 生 
态 环境 保护 意识 ,促进 区 域 多 种 生态 系统 服务 良性 
发 展 。 


3 讨论 
3.1 生态 系统 服务 空间 差异 性 分 析 
生态 系统 过 程 是 自然 和 人 类 活动 共同 作用 的 


结果 ,如 何 构建 合理 的 模型 来 模拟 区 域 生态 系统 服 
务 是 深刻 了 解 生态 过 程 的 关键 ,本 文通 过 不 同 的 生 
态 系 统 服 务 评估 模型 得 到 4 种 生态 系统 服务 的 空间 


(a) NPP 优 化 区 


(b) 产 水 量 优化 区 


格局 。 为 了 探究 研究 结果 的 准确 性 ,将 本 文 的 研究 
结果 和 其 他 学 者 进行 对 比 。 刘 宇 等 外 通 过 用 Pen- 
man- Monteith 公式 计算 泾 河流 域 近 70 a 的 径流 变 
化 ,以 及 黄 晨 璐 等 “在 统计 数据 的 基础 上 ,运用 相 
关 的 模型 来 模拟 泾 河流 域 径流 状况 ,在 时 空 尺 度 上 
和 本 研究 结果 具有 相似 性 ;在 NPP 和 土壤 保持 服务 
的 时 空 分 布 方面 与 卞 鸿雁 等 “以 及 Zheng 等 ”的 结 
果 具 有 相似 性 ,都 在 年 际 间 趋 于 稳定 上 升 的 趋势 ， 
并 且 时 空 差异 显著 。 此 外 ,根据 本 文 研究 结果 ,4 种 
生态 系统 服务 的 优化 区 域 都 集中 在 黄土 高 原 和 泾 
河流 域 北部 的 山区 ,制约 区 域 发 展 的 关键 因子 主要 
是 土壤 .植被 和 降水 ,而 关键 因子 对 于 评估 区 域 4 种 
生态 系统 服务 有 着 至 关 重 要 的 影响 。 泾 河流 域 产 
水 服务 的 空间 变化 特征 和 土地 利用 、 蒸 散发 和 植被 
水 分 利用 系数 相关 ,2000 一 2020 年 区 域 耕 地 面积 减 
少 , 和 森林 和 草地 面积 增加 ,进而 增加 区 域 植被 水 分 


(c) 农业 生产 力 优化 区 


图 例 

O 县 区 边界 
生态 系统 服务 优化 区 
NPP 优 化 区 

国 | 产 水 量 优 化 区 

辐 ] 农业 生产 力 优 化 区 
mm 土壤 保持 优化 区 
综合 优化 区 


0 50 km 
L_ 


图 5 综合 生态 系统 服务 优化 区 


Fig. 5 Integrated ecosystem service optimization zone 
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利用 系数 ,植被 类 型 从 禾 本 科 植 物 向 森林 植被 转化 
所 引起 。2020 年 相 较 于 2000 年 ,降水 量 增加 ,而 部 
分 区 域 径流 量 减少 ,暗示 区 域内 潜在 蔡 散 发 则 逐渐 
增 大 ,表明 区 域 径流 量 对 于 降水 和 潜在 蒸 散发 较为 
敏感 于。 产 水 和 降水 的 空间 分 布 较为 一 致 ,主要 取 
决 于 气候 和 地 形 因素 相关 ,西北 部 属于 黄土 高 原 


各 种 方案 下 的 服务 变化 ,并 提供 优化 建议 “"”。 然 
而 ,土地 利用 并 不 是 唯一 影响 生态 系统 服务 的 因 
素 ,而 各 因素 对 不 同 生态 系统 服务 的 影响 程度 也 是 
不 同 的 。 本 文 在 贝 叶 斯 网 络 的 基础 上 ,提出 了 关键 
子 集 法 ,对 泾 河流 域 生 态 系统 服务 的 空间 格局 进行 
了 优化 ,研究 结果 为 研究 区 生态 环境 保护 区 划 提 供 


区 , 受 副热带 高 气压 控制 ,干旱 少雨 , 产 水 量 也 较 
少 。 泾 河流 域 大 多 区 域 经 济 相对 落后 ,居民 的 生活 
需求 对 区 域 生态 系统 服务 依赖 性 很 强 ,农业 生产 力 
较 强 的 区 域 ,居民 能 够 满足 生活 ,少数 区 域 粮 食 供 
给 能 力 不 足 ”。 北 部 区 域 以 草地 为 主 , 牧 业 长 期 集 
中 与 此 ,长 期 的 粗放 经 营 ,会 造成 土地 严重 的 沙化 ， 
进而 影响 区 域 的 农业 生产 力 以 及 NPP 服 务 ,造成 两 
种 生态 系统 服务 在 区 域 上 呈现 南北 差异 化 明显 。 
而 随 着 退耕 还 林 封山育林 .飞播 .人 工 植树 等 生态 
措施 的 实施 ,促使 区 域内 林地 面积 增加 ,提升 区 域 
的 NPP 和 水 土 保持 能 力 。 本 人 研究 建议 对 生态 系统 
服务 优化 区 实施 提升 肉 类 的 产量 以 及 薪 荣 供给 的 
手段 等 生态 措施 和 生态 补偿 政策 ”” ,不仅 可 以 提 
升 居民 的 收入 和 生活 质量 ,而 且 还 可 以 保障 原 有 的 
土壤 保持 和 NPP 等 生态 系统 服务 的 增强 ,促进 区 域 
生态 经 济 和 谐 发 展 。 为 了 深入 人 研究 生态 系统 服务 
的 区 域 差 异性 ,不 仅 要 定性 分 析 , 还 要 利用 模型 进 
行 定量 的 研究 ,未 来 应 该 通过 多 视角 、 多 指标 对 生 
态 系统 服务 驱动 因素 进行 定量 分 析 ,促使 研究 结 
更 加 精准 。 
3.2 生态 系统 服务 空间 格局 优化 

贝 叶 斯 网 络 可 以 将 多 元 知识 进行 图 解 可 视 化 
的 概率 模型 , 它 不 仅 包括 各 个 网 络 节 点 之 间 的 因果 
相关 性 ,而 且 结 构 还 比较 灵活 , 既 可 以 优化 现在 的 
空间 格局 ,还 可 以 预测 未 来 的 空间 格局 分 布 。 李 婷 
等 "通过 贝 叶 斯 网 络 和 固 碳 服务 相 结合 ,来 探究 固 
碳 服务 的 格局 优化 问题 ; 兽 莉 等 ”通过 不 同 发 展 情 
景 水 源 涵养 的 状态 概率 分 布 ,最 终 给 出 水 源 涵养 服 
务 的 空间 格局 优化 策略 。 但 是 贝 叶 斯 网 络 是 具有 
不 确定 性 的 ,主要 包括 结构 不 确定 性 、 输 入 数据 不 
确定 性 和 参数 不 确定 性 ” 。 本 文 将 产 水 量 .NPP 、 土 
壤 保 持 和 农业 生产 4 种 生态 系统 服务 和 贝 叶 斯 网 络 
相 结 合 , 探 讨 了 不 同 生态 系统 服务 以 及 综合 优化 区 
Jk ,为 优化 区 的 保护 和 生态 经 济 和 谐 发 展 黄 定 基 
础 。 目 前 ,对 生态 系统 服务 空间 格局 优化 的 研究 主 
要 集中 在 围绕 土地 利用 模式 的 变化 设计 方案 ,预测 


决策 支持 ,丰富 生态 系统 服务 空间 格局 优化 的 研究 
方法 。 但 是 自然 和 人 类 活动 共同 影响 着 区 域 的 生 
态 环 境 ,本文 仅 考虑 降水 量 气温 .土地 利用 等 因素 
对 生态 系统 的 影响 是 不 够 的 ,未 来 应 该 结合 更 多 的 
影响 指标 ,比如 人 类 活动 .生态 政策 以 及 社会 经 济 
等 因素 来 构建 多 情景 的 模式 ,将 情景 分 析 和 贝 叶 斯 
网 络 相 结合 ,更 加 全 面 的 为 区 域 空间 格局 优化 提供 
依据 。 


4 结论 


(1) SWAT 模 型 在 泾 河流 域 有 较 好 的 适用 性 。 
校准 后 的 SWAT 模 型 对 径流 量 的 模拟 结果 与 实测 值 
相 比 ,具有 较 高 的 决定 性 系数 (RR>0.6) .纳什 效 率 系 
数 (NSE>0.5) ,表明 该 模型 能 较为 准确 地 模拟 流域 
生态 水 文 过 程 ,为 进一步 评估 产 水 服务 提供 保障 。 

(2) 2000 一 2020 年 泾 河流 域 4 种 生态 系统 服务 
的 时 空 差 异性 较为 显著 。 在 时 间 尺 度 上 , 随 着 退耕 
还 林 还 草 以 及 居民 生态 环保 意识 的 提升 ,促使 4 种 
生态 系统 服务 都 呈现 波动 中 上 升 的 趋势 ;而 在 空间 
尺度 上 ,NPP 农业 生产 力 以 及 土壤 保持 服务 变化 较 
小 ,呈现 较为 稳定 的 趋势 ,而 产 水 量 和 区 域 气候 有 
着 直接 的 关系 ,导致 产 水 服务 在 区 域内 呈现 显著 
变化 。 

(3) 贝 叶 斯 网 络 模型 是 一 种 强大 的 神经 网 络 工 
具 , 可 以 很 好 地 与 生态 系统 服务 相 结合 ,在 生态 系 
统 服务 空间 优化 过 程 中 发 挥 重要 作用 。 本 研究 通 
过 对 4 种 生态 系统 服务 优化 区 域 进行 到 加 分 析 ,发 
现 综合 优化 区 域 集中 在 彭 阳 县 的 中 部 和 西南 部 以 
及 环 县 的 零星 区 域 ,对 区 域 应 该 采取 相应 的 生态 环 
境 保护 措施 ,要 坚持 保护 与 建设 并 举 ,保护 优先 的 
原则 ,进而 促进 区 域 多 种 生态 系统 服务 良性 发 展 。 
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Spatial pattern optimization of ecosystem services based on Bayesian networks: 
A case of the Jing River Basin 


YU Yuyang", LI Jing, ZHOU Zixiang, TANG Chengyan' 
(1. School of Geography and Tourism, Shaanxi Normal University, Xi "an 710119, Shaanxi, China; 2. School of Tourism, Henan 
Normal University, Xinxiang 453007, Henan, China; 3. College of Geomatics, Xi an University of Science and Technology, 
Xi'an 710054, Shaanxi, China) 


Abstract: Changes in regional land use patterns due to rapid urbanization affect the quality of the ecological 
environment and the spatial pattern of ecosystem services. Therefore, optimizing ecosystem services is crucial. 
Using the spatial assessment of net primary production, crop production, soil conservation, and water yield 
services in the Jing River Basin, northwest China from 2000 to 2020 as the basis, this study combined Bayesian 
networks and ecosystem services and employed a subset of key variables and a visualized subset of optimal states 
to determine the areas necessitating the optimization of the aforementioned ecosystem services for regional 
economic and ecological harmonious development. Results showed that (1) the SWAT model can accurately 
simulate the regional runoff and generate a high coefficient of determination (R>0.6) and a Nash-Sutcliffe efficiency 
coefficient (NSE>0.5) when comparing simulated and observed values, thus providing a guarantee for the further 
assessment of water yield services. (2) The spatial and temporal variations of the four ecosystem services in the 
Jing River Basin from 2000 to 2020 were substantial. All four ecosystem services showed an upward fluctuation 
trend on the temporal scale and a stable trend on the spatial scale. (3) Overlay analysis of the optimized areas for 
the four ecosystem services revealed that the integrated optimized areas were concentrated in the central and 
southwestern parts of Pengyang County, Ningxia Province and the scattered areas of Huan County, Gansu 
Province. This study is important for guiding region optimization for the sustainable management and reducing 
the degradation of ecosystems. 


Key words: ecosystem services; Bayesian networks; spatial and temporal pattern optimization; Jing River Basin 


